An Electromagnetic GL Double Layered Cloak 
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In this paper, we propose a new electromagnetic (EM) GL double layered cloak. The GL double 
layered cloak is consist of two sphere annular layers, Ri < r < R2 and R2 < r < R3. Two type 
cloak materials are proposed and installed in the each layer, respectively. The outer layer cloak of 
the GL double layered cloak has the invisible function, the inner layer cloak has fully absorption 
function. The GL double layered metamaterials are weak degenerative and weak dispersive. When 
the source is located outside of the GL double layered cloak, the excited EM wave field propagation 
outside of the double layered cloak is as same as in free space and never be disturbed by the cloak; 
also, the exterior EM wave can not penetrate into the inner layer and concealment. When local 
sources are located inside of the GL double cloaked concealment with the normal EM materials, 
the excited EM wave is propagating under Maxwell equation governing, it is complete absorbed 
by the inner layer cloak of GL double cloak and never propagate to outside of the inner layer of 
the GL cloak, moreover, the EM wavefield in concealment never be disturbed by the cloak. The 
GL doubled layered cloak is a robust cloak and has complete and sufficient invisibility functions. 
Its concealment is the normal electromagnetic environment. Our EM GL double layered cloak is 
different from conventional common cloak. The 3D GL EM modeling simulations for the double 
layered cloak are presented. The GL method is an effective physical simulation method and is 
fully different from the conventional methods. It has double abilities of the theoretical analysis 
and numerical simulations to study the cloak metamaterials and wide materials and field scattering 
problem in physical sciences. 

PACS numbers: 13.40.-f, 41.20.-q, 41.20.jb,42.25.Bs 



I. INTRODUCTION 



Using the 3D GL EM modeling [1-2] and GL inversion 
[3] , we propose an EM double layered cloak in this paper 
which is called as GL double layered cloak. The single 
layer cloak proposed by Pendry et al. [4] is named as 
PS cloak. The GL double layered cloak is consist of two 
sphere annular layers, R\ < r < R2 and R2 < r < R$. 
Two different type cloak materials are proposed and in- 
stalled in the each layer, respectively. The outer layer 
cloak of the GL double layered cloak has the invisibility 
function, while the inner layer cloak has fully absorp- 
tion function to absorb the EM wave excited from local 
sources inside of the concealment. When the source is 
located outside of the GL double layered cloak, the ex- 
cited EM wave field propagation outside of the double 
layered cloak is as same as in free space and never be dis- 
turbed by the cloak; also, the exterior EM wave can not 
penetrate into the inner layer and concealment. When 
local sources are located inside of the GL double cloaked 
concealment with the normal EM materials, the excited 
EM wave is propagating under Maxwell equation govern- 
ing, it is complete absorbed by the inner layer cloak of 
GL double cloak and never propagate to outside of the 
inner layer of the GL cloak, moreover, the EM wavefield 
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in concealment never be disturbed by the cloak. 

The GL doubled layered cloak is a robust cloak and has 
complete and sufficient invisibility functions. Its conceal- 
ment is the normal electromagnetic environment. Our 
EM GL double layered cloak is different from conven- 
tional common cloak. The GL double layered metama- 
terials are weak degenerative and dispersive. 

A double layer cloth phenomenon to prevent the 
GILD inversion [5-7] detection has been observed 

in pa per [8] in 2001 which is published in SEG online; 

http://www.segdl.org/journals/doc/SEGLIB-home/dci/searchDCL. 
We developed a novel and effective Global and Local 
field (GL) modeling and inversion to study the metama- 
terials, periodic photonic crystals and condense physics 
etc. wide physical sciences. 3D GL EM modeling and 
inversion [9] and computational mirage [10] have been 
presented in PIERS 2005 and published in proceeding of 
PIERS 2005 in Hangzhou, which can be downloaded from 
http://piers.mit.edu/piersproceedings/piers2k5Proc.php. 
The analytical method and numerical method for phys- 
ical sciences are developed separately in history. The 
GL method consistent combines the analytical and 
numerical methods together. The GL method does not 
need to solve large matrix equation, it only needs to 
solve 3x3 and 6x6 matrix equations. Moreover, the 
GL method does not need any artificial boundary and 
not need PML absorption condition on it to truncate 
the infinite domain. The Finite Element method (FEM) 
and Finite Difference (FD) method have numerical 
dispersions which confuse and contaminate the physical 
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dispersion from the interaction between the field and 
materials. The frequency limitation is a difficulty of 
FEM and FD method. 

The GL method is a significant scattering process 
which reduces the numerical dispersion and is suitable to 
simulate physical wavefield scattering in the materials, in 
particular, for dispersive materials. Born Approximation 
is a conventional method in the quantum mechanics and 
solid physics, however, it is one iteration only in whole 
domain which is not accurate in the high frequency and 
high contrast materials. The GL method divides the do- 
main as a set of small sub domains or sub lattices. The 
Global field is updated by the local field from the in- 
teraction between the global field and local subdomain 
materials successively. Once all subdomain materials are 
scattered, the GL field solution is obtained which is much 
more accurate than the Born approximation. 

Moreover, the GL method can be meshless, including 
arbitrary geometry subdomains, such as rectangle, cylin- 
drical and spherical coordinate mixed coupled together. 
It is full parallel algorithm. These advantages of the GL 
method to overcome historical difficulties have been de- 
tailed described in the paper [I]. The theoretical foun- 
dation of the GL method is described in the paper [2] 

We have used the 3D GL modeling [1-2] and inversion 
[3] to simulate many cloak metmaterials, nanometer ma- 
terials, periodic photonic crystals etc. When the point 
source is located outside or inside of the various geome- 
try cloaks, the 3D GL EM modeling simulations for the 
EM wavefield propagation through the cloaks have been 
done. These simulations show that the GL method is fast 
and accurate. We have submitted a paper titled the 3D 
GL EM modeling to simulate single layer cloaks to PRE 
[11]. 

In this paper, the 3D GL EM modeling simulations of 
the EM wave field propagation through the new GL dou- 
ble layered cloak is presented. When the local sources 
are located outside or inside of the outer layer cloak, EM 
wave propagation through outer layer cloak and never 
penetrate into the inner layer and the concealment, i.e. 
r < i?2- The exterior EM wavefield propagation outside 
of GL double layered cloak never be disturbed by the 
cloak. The outer layer cloak has the invisibility function. 
When the local sources inside of the GL double layer 
cloaked concealment with normal materials, the excited 
EM wave normally propagating under Maxwell equation 
governing, the EM field is complete absorbed by the in- 
ner layer cloak and can not propagate outside of the inner 
cloak. The EM environment in the GL double cloaked 
concealment is normal, in which there exist Maxwell EM 
wave field excited by nonzero local sources, have no re- 
flection from the boundary r = R\, and never propagate 
outside of boundary r = R2 

By using the 3D GL EM modeling [1-2] simulation 
and its theoretical analysis, we found and verified a phe- 
nomenon that there exists no Maxwell EM wave field can 
be excited by nonzero local sources inside of the single 
layer cloaked concealment with normal materials. Our 



GL double layered cloak overcomes the drawback and 
difficulty in the single cloak. Pcndry et al. in paper 
4] used a coordinate transformation and ray tracing to 
propose the annular cloak in which the ray being bending 
and re direction around central sphere object and can not 
penetrate into it. The cloak device like empty and does 
not disturb exterior wave field. There are several other 
papers to simulate the exterior plane wave propagation 
through the cloak [12-14]. Cummer et al. in paper [l2j 
proposed numerical simulations by using the COMSOL 
Multiphysics finite element-based electromagnetics solver 
for the 2D plane wave propagation through cylindrical 
cloak. Chen et al. proposed the Mie analytical TEM 
model to simulate the plane wave through the spherical 
cloak [TI|. Argyropoulos et al. proposed a dispersive fi- 
nite difference method in time domain (FTFD) in [14] 
to simulate 2D TEM plane wave field through cylindri- 
cal cloak, in which authors considered the difficulty of 
conventional FDTD scheme for dispersive materials. In 
papers [l^] and [l4| , authors introduced many papers for 
cloak research works. Because the plane wave is excited 
by plane source which can not be located inside of the 
cloak or concealment. To study the EM wave excited 
from local sources inside of the cloaked concealment is 
absent from these papers. In paper [15] and [16], authors 
studied the effect on invisibility of active devices inside 
the cloaked region. Author in [16] stated that "when 
these conditions are overdermined, finite energy solutions 
typically do not exist." 

We use 3D GL method to do many simulations for 
studying the behavior of EM field excited inside of the 
single layer cloaked concealment. These simulation are 
divergent or become chaos when the EM wave propa- 
gates to the inner boundary of the single layer cloak. 
Our statement is that "there exists no Maxwell EM wave 
field can be excited by nonzero local sources inside of 
the single layer cloaked concealment with normal mate- 
rials". The detailed proof and 3D GL simulations are 
presented in this paper. Before the practice production 
of the single layered cloak, the electromagnetic field en- 
vironment inside of the concealment can not be studied 
in physical experiment. Our GL double layered cloak 
proposed in this paper overcomes the drawback and dif- 
ficulty of the single layer cloak and avoid the disputing on 
the EM phenomenon inside of the single layered cloaked 
concealment with normal materials. The GL double lay- 
ered EM cloak metamaterials inventive and fabrication 
technology right and 3D GL EM modeling software are 
patented by GL Geophysical Laboratory. We thanks to 
GL Geophysical Laboratory to approve us to publish the 
GL modeling method, GL double layered cloak theory, 
and simulations. 

We describe this paper in the following order: The in- 
troduction is described in Section 1. In Section 2, we 
propose a GL double layered cloak materials. The EM 
integral equations are presented in Section 3. The 3D GL 
EM modeling are described in Section 4. The theoreti- 
cal analysis of properties and functions of the GL double 
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layered cloak are proposed in Section 5. The simulations 
of the EM wave propagation through the GL double lay- 
ered Cloak by using the GL EM modeling are presented 
in Section 6. The advantages of the GL double layered 
cloak is presented in Section 7. In Section 8, we conclude 
this paper. 



The GL inner layer cloak material [D] GLI = diag [ei,pi] 
in (1) on the Qgli an d GL outer layer cloak material 
[D] GLO = diag [e , p, ] in (2) on Qglo are assembled into 
the GL anisotropic double layered cloak material on the 
domain Qgl- The domain SIql with the metamaterial 
[D] GL in (4) is called as the GL double layered cloak. 



II. GL DOUBLE LAYERED CLOAK 
MATERIALS 



III. 3D ELECTROMAGNETIC INTEGRAL 
EQUATION 



A. GL Inner Layered Cloak Anisotropic Material 

On the inner sphere annular domain, Qgli — 
{r : Ri < r < R 2 } , by the GL EM modeling and inver- 
sion [1-3], we propose an anisotropic material as follows, 

[ d \gli = dia 9 , 
e» = diag \e r ,%, ee,», £#,»] £o, 
fii = diag [/ir,», M0,*> M<A,i] Mo, 

^2 ) 



R 2 -r 2 



£0 



— -^2 

R 2 -r 2 R 2 -r 2 ' 



The ftcLi is called as GL inner layered cloak, the mate- 
rials, [D]gli = diag [Si, fii] in (1), are the anisotropic GL 
inner layered cloak material tensor. 



B. GL Outer Layered Cloak Anisotropic Material 

Let the outer sphere annular domain Qglo — 
{r : R 2 < r < R 3 } be the GL outer layered cloak with 
the following anisotropic GL outer layered cloak materi- 
als, 



[°\glo = dia 9 [e ,Mo] > 
e = diag [e r ,o, £e,o, £</>,o] £o, 
fio = diag [/i rjD , fj, s>0 , fi^ ] fj, , 



- _,, _ R 3 r'-Rl 

e r ,o - Mr,o - r —pr- ^Zri 
£e.o = Me.o = £cf>,o = ^4>.o 

_ R 3 r 



(2) 



The GL inner cloak Qqli domain and GL outer cloak 
Qglo domain are bordering on the interface surface 
r = R 2 . We assemble the Qgli as the inner layer sphere 
annular domain and Qglo as the outer layer sphere an- 
nular domain and make them coupling on their interface 
boundary annular surface r — R 2 as follows, 

&GL = &GLI U ^GLO 

= {r:R 1 <r<R 2 }{J{r:R 2 <r<R 3 } (3) 
= {r : Ri < r < R 3 } , 

and offer the coupled anisotropic dielectric and suscepti- 
bility tensor [D] GL on the Qgl as follows, 



[D] GL = { 



[D] 

GLI ' r ^ ^GLI 
[D] GLO ' r e ^GLO- 



(4) 



The 3D EM integral equation in frequency domain has 
been proposed in authors' papers [1] and [2]. In this 
paper, we proposed the EM integral equation in time 
domain as follows: 



(5) 
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+ fGi M H 
(1) a ' 


(r',r 


,t) * t S[D(r')} 


' E b {r',t) ' 
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' E(r,t) ' 
H(r, t) 




' E b (r,t) ' 
_H b (r,t) 









n 



E(r',t) 
H(r',t) 



dr'. 



In the EM integral equation (5), 



G J d M H (r',r,t) 



E J {r',r,t) H J {r',r,t) 
E M (r' , r, t) H M {r',r,t) 



(6) 



(7) 



and b is the EM Green's tensor in the background 

medium, where, E(r, t) is the electric field, H(r, t) is the 
magnetic field, E b (r, t) and H b {r, t) is the incident electric 
and magnetic field in the background medium, E J (r' , r, t) 
is electric Green's tensor, H J (r', r, t) is magnetic Green's 
tensor, they are excited by the point impulse current 
source, E M (r' , r, t) and H M (r' , r, t) are electric and mag- 
netic Green's tensor, respectively, they are excited by the 
point impulse magnetic moment source, *t is convolution 
with respect to t, S [D] is the electromagnetic material 
parameter variation matrix, 



S[D] 





(a(r) 



6D 22 = (fi(r) - » b I)f 



£bI)-§- t , 



(8) 





5D 22 

a b I) + (e(r) 
i df 

SDn and SD 22 are a 3 x 3 symmetry, inhomogeneous di- 
agonal matrix for the isotropic material, for anisotropic 
material, they are an inhomogeneous diagonal or full ma- 
trix, / is a 3 x 3 unit matrix, <r(r) is the conductivity 
tensor, e(r) is the dielectric tensor, p,(r) is susceptibility 
tensor which can be dispersive parameters depend on the 
angular frequency w, a b is the conductivity, e b is the per- 
mittivity, /i b is the permeability in the background free 
space, £1 is the finite domain in which the parameter vari- 
ation matrix 6 [D] ^ 0, the (e(r) — s b I)E is the electric 
polarization, and (fi(r) — fi b I)H is the magnetization. 
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IV. 3D GL EM MODELING 

We propose the GL EM modeling based on the EM 
integral equations (1) and (2) in the time space domain. 

(3.1) The domain fi is divided into a set of N sub 

N 

domains, {fi^}, such that fi = (J fi^. The division can 

fe=i 

be mesh or meshless. 

(3.2) When k = 0, let E (r,t) and H (r,t) are the an- 
alytical global field, E^(r', r, t), H£{r', r, t), E^(r', r, t), 
and Hq 1 (r' ,r,t) are the analytical global Green's ten- 
sor in the background medium. By induction, suppose 
that £ fc _i(r,t), H k -i(r,t), E J k _ x {r> >,i), H^r' ,r,t), 
-^feli( r ') r )*)i ancl H^ 1 (r / ,r,t) are calculated in the 
(k — l) th step in the subdomain fifc-i- 

(3.3) In {fife}, upon substituting E k -i(r,t), 
ff fc _i(r,t), E(_ l {r J ,r,t), H^r' ,r,t), E^r' ,r,t), 
and H^_ 1 (j J ,r,t) into the integral equation (1), the EM 
Green's tensor integral equation (1) in fife is reduced 
into 6x6 matrix equations. By solving the 6x6 matrix 
equations, we obtain the Green's tensor field E^(r', r, t), 
Hi (r' ,r,t), Ejf (r' , r, t) , and fff (r' , r, t) . 

(3.4) According to the integral equation (1), the elec- 
tromagnetic field Ek(r,t) and H k (r,t) are updated by 
the interaction scattering field between the Green's ten- 
sor and local polarization and magnetization in the sub- 
domain fife as follows, 



(9) 



" E k {r,t) ' 




' B*-i(r,t) " 


_H k (r,t)_ 




_ (r,t) _ 



+ / 



E J k (r>,r,t) m(r',r,t) 



E^(r',r,t) Hjf(/,r,t) 



t8[D(f)] 



E k -i(r',t) 



dr> 



(3.5) The steps (3.2) and (3.4) form a finite iteration, 
k = 1, 2, • • • , N, the E N (r, t) and H N (r, t) are the elec- 
tromagnetic field of the GL modeling method. The GL 
electromagnetic field modeling in the time space domain 
is short named as GLT method. 

The GL EM modeling in the space frequency domain 
is proposed in the paper 0, we call the GL modeling in 
frequency domain as GLF method. 



V. THEORETICAL ANALYSIS OF 
INTERACTION OF THE EM WAVE FIELD 
THROUGH THE CLOAKS 

A. Theoretical Analysis Of Interaction Of The EM 
Wave Field Through The GL Double Layered Cloaks 

We propose the theoretical analysis of the interaction 
between the EM wave and GL cloaks in this section. 

Statement 1: Let domain ^Iql m (3) and the meta- 
material Dql in (4) be GL double layered cloak, and 
s=Eh,fJ- — fJ-b be basic permittivity and permeability, 



respectively, inside of the central sphere concealment 
\r'\ < Ri and outside of the GL cloak \r'\ > i?3, we 
have the following statements: (1) provide the source is 
located inside of the concealment of GL double layered 
cloak, | fa | < Ri, the excited EM wave field propagates 
inside of the concealment and never be disturbed by the 
cloak; (2) provide the source is located inside of conceal- 
ment or inside of the inner layer of the GL double layered 
cloak, \r s \ < R2, the EM wave field is vanished outside 
of the inner layer of GL cloak and is always propagating 
to the boundary r = R2 and is absorbed by the bound- 
ary r = i?2- (3) provide the source is located outside of 
the GL double layered cloak, |r s | > i? 3 , the excited EM 
wave field outside of the double layered cloak is as same 
as in free space and never be disturbed by the double 
layered cloak; (4) provide the source is located outside of 
double layered cloak or located inside of the outer layer 
of GL cloak, \f s \ > R2, the excited EM wave field never 
propagate into the inner layer of GL cloak and the con- 
cealment. 



B. There Exists No Maxwell EM Wavefield Can 
Be Excited By Nonzero Local Sources Inside Of The 
Single Layered Cloaked Concealment With Normal 
Materials 



Statement 2: Suppose that a 3D anisotropic inhomo- 
geneous single layered cloak domain separates the whole 
3D space into three sub domains, one is the single lay- 
ered cloak domain Q c ik with the cloak material; the sec- 
ond one is the cloaked concealment domain fi COTl ; with 
normal EM materials; other one is the free space outside 
of the cloak. If the Maxwell EM wavefield excited by a 
point source or local sources outside of the concealment 
ficon; is vanished inside of the concealment fi corl j, then 
there is no Maxwell EM wave field excited by the local 
sources inside of the cloaked concealment il con i ■ 

The Maxwell EM wavefield is the EM wave field which 
satisfies the Maxwell equation and tangential continuous 
interface boundary conditions. We call the Maxwell EM 
wavefield as the EM wavefield and use inverse process to 
prove the statement 2 as follows: Suppose that there ex- 
ists Maxwell EM wavefield excited by the local sources 
inside the concealment with the normal materials, the 
wavefield satisfies the Maxwell equation in the 3D whole 
space i? 3 which includes the anisotropic inhomogeneous 
cloak domain fi c ;fe and concealment fi con /, and satisfies 
the tangential continuous interface conditions on the in- 
terface boundary surface Si and S2. The Si is the in- 
terface boundary surface between the cloak domain fi c ;fe 
and the concealment fi cori ;, it also is the inner boundary 
surface of the cloak domain fi c /fe . The S2 is the interface 
boundary surface between the cloak domain fi c ;fe and the 
free space, it also is the outer boundary surface of the 
cloak domain fi c ;fe. 

Let R c — R 3 — fi cZfe (J fi con ;, R d = R 3 — fi co „ ; , and by 
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the EM integral equation (1), the EM wave field satisfies conditions of G^j^ (r',r,t), the term in right hand side 



~ E(r,t)~ 




' E b (r,t) ' 


_H(r,t)_ 




H b (r,t)_ 



E b (r',t) 
H b (r',t) 



dr', 



(10) 



/ Gi M H (r',r,t) H S[D] 

where G^^ (r' , r, i) is the EM Green's tensor, its compo- 
nents E J , H J , E M , and H M (r' , r, t) arc the EM Green's 
function on fl c ik{jQ C oni{J R c , excited by the point im- 
pulse sources outside of the concealment, r G Rd- By the 
assumptions, G^^(r',r, t) exists on fl c ik \J £l CO ni U R c 

and when r' G Sl C oni, G~e% ( r ' > r > *) = 0- The integral 
equation (10) becomes to 



' E(r,t)' 




' E b (r,t) ' 


_H(r,t)_ 




_H b {r,t) 



+ I G J E M H (r>,r,t)* t 6[D] 



E b (r',t) 
H b (r',t) 



dr'. 



(11) 



We consider the Maxwell equation in Rd, the virtual 
source is located r, r G Rd and the point source is located 
r s , r s G fi co „ ; and r s R d , we have 



Vx 



-Vx 



[D] G J i M H (r',r,t)+I8(r',r)S(t), 



(12) 



and 



-Vx 
= [D b ] 



Vx 

E b ' 
H b 



E h 
H b 



(r',r s ,t) 



(r',r s ,t) 



(13) 



By using [E b (r, t) , H b (r, i)] to convolute (12), and 
G E ^ (r' ,r,t) to convolute (13), to subtract the second 
result equation from the first result equation and make 
their integral in Rd, and use integral by part and make 
some manipulations, we can prove 



E b (r,t) 
H b (r,t) 

+ I G J E %{r',r,t) H 5[D] 



E b {r',t) 
H b (r',t) 



dr' 



(14) 



§ G^(r',r,i)< 

Si 



E b {r',t) 
H b (r',t) 



dS, 



where (S>j is cross convolution. From the assump- 
tion of the statement 2 that "the Maxwell EM wave- 
field excited by a point source or local sources out- 
side of the concealment VL CO ni is vanished in inside of 
the concealment f2 con z", and virtual source r is located 
outside of the concealment, r G Rd, if r' G Q CO ni, 
Ge^h ( r 'i r it) = 0- By continuity, when r' G Si, we have 
G E ^(r',r,t) = 0. By tangential continuous interface 



of (14) is vanished, we have 



E b (r,t) 
H b {r,t) 

f n J,M 



+ 



(r',r,t) * t 5[D] 



E{r',t) 
H(r',t) 



dr' = 0. 



(15) 



Upon substituting integral equation (15) into the inte- 
gral equation (11), we have 



E(r,t) 
H(r,t) 



= 0. 



(16) 



From the continuous property of the EM wave field, 
we obtain the following over vanish boundary condition 
on the boundary Si of the concealment Q CO ni, we have 



E(r, r s ,t) 
H(r,r s ,t) 



= 0, 



(17) 



Si 



where r s denotes point source location inside of the con- 
cealment, r is EM field receiver point, r G Si. Because 
the EM wave field is excited by local sources inside of 
the concealment domain £l con i, it satisfies the following 
Maxwell equation, 



Vx 



-Vx 





' E ' 




H 



E 
H 



(r',r s ,t) 
(r',r s ,t) + Q(r',r s ,t). 



(18) 



where [D con i] = diag [e r e u r n o] (d/dt) with the nor- 
mal EM material parameters, e r > 1 and \i r > 1 are 
relative EM parameters, £o is basic permittivity and /xo 
is basic permeability, r s G Qconi 

is the local source loca- 
tion, Q(r' , r s ,t) is the nonzero local source inside of Q con i- 
Let G E ^j conl (r' ,r,t) be Greens tensor which satisfies 



Vx 
-Vx 



G 



J,M 

E,H.conl 



(r',r,t) 



[D con i]G^ H nl {r' >, t) 



(19) 



G 



+IS(r',r)5(t) 

By using [E (r, t) , H (r, t)} to convolute (19), and 
(r',r,t) to convolute (18), to subtract the sec- 



J,M 

E,H.conl 



ond result equation from the first result equation and 
make their integral in O co „;, and use integral by part 
and make some manipulations, we have 



E{r,r s ,t) 
H(r, r s ,t) 
I G 



+ § 

dn cor 



G 



J,M 

E,H,conl 
J,M 



(r',r, t) * t Q(r',r s ,t)dr' 



(20) 



E,H,con, 



i(r',r, t) ® t 



E(r',r s ,t) 
H(r',r s ,t) 



dr', 



<8>t denotes the cross convolution, and d£l c 
cause the over vanished condition (17), 



Si. Be- 



E(r, r s ,t) 
H(r, r s ,t) 



= 0, 
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we have 



E(r, r s 
H(r, r s 



t) 
t) 

E,H-< 



t {r',r,t) * t Q(r',r s ,t)dr'. 



(21) 



Because G 



J,M 

E ,H \conl 



(r',r,t) ^ and Q(r',r s ,t) ^ 0, so, 



E(r,r s ,t) 
H(r, r s ,t) 



From the continuity of the EM 
EM wave field (22) results that 



^0. 



(22) 



wavefield, the nonzero 



E(r, r s ,t) 
H(r,r s ,t) 



7^0. 



(23) 



The 
ary 
the 
ary 



EM wavefield is nonzero on the bound- 
s'! in (23) is an obvious contradiction with 
same EM wavefield is zero on the bound- 
5*i in (17). Therefore, we proved that 



VI. THE GL EM MODELING SIMULATIONS 
OF THE EM WAVE FIELD THROUGH THE GL 
DOUBLE CLOAKS 



A. The Simulation Model of The GL Double 
Layered Cloak 

The simulation model: the 3D domain is 
[—0.5m, 0.5m] x [—0.5m, 0.5m] x [—0.5m, 0.5m], the 
mesh number is 201 x 201 x 201, the mesh size is 0.005m. 
The electric current point source is defined as 



5{r - r s )S(t)e, 



(29) 



where the r a denotes the location of the point source, the 
unit vector e is the polarization direction, the time step 
dt = 0.3333 x 10 -10 second, the frequency band is from 
0.05 GHz to 15 GHz, the largest frequency / = 15GHz, 
the shortest wave length is 0.02m. The EM GL double 
layered cloak Qql = Qgli U ^glo is consist of the dou- 



there exists no Maxwell EM wave field can be excitedjlg^r^ical annular n GLI and VL G lo with the center in 
nonzero local sources inside of the single layered cloa4fed or igi n an d interior radius R x = 0.22m, meddle ra- 



concealment. For more simplicity to derive the nonzero 
EM wavefield (22) from the integral expression of the 
EM wave field (21), let the source is point impulse 
current source with polarization direction in x, i.e., 



Q(r,r S) t)=6(r-r s )6(t)x. 



(24) 



dius i?2 = 0.3m. and exterior radius i?3 — 0.35771. The 
cloak is divided into 90 x 180 x 90 cells. The spherical 
coordinate is used in the sphere r < i?3, the Cartesian 
rectangular coordinate is used in outside Qgl to mesh 
the domain. 



Upon substituting (24) and e r = 1.0 and \x r — 1.0 into 
the (21), we have 



E(r, r„t) 
H(r, r„t) 



E J x {r,r s ,t) 



E£(r,r„t) 



E xx (r, r s ,t) 
E xy (r, r s ,t) 
E xz (r, r s ,t) 



E xx (r, r s ,t) 

1 d 2 S(t-^/ejI\r-r B \) 

Stt 2 £ dx' 2 \r— r s \ 

! Qi &{t-^i\r-r,\) 



\t—T s I 



(25) 



(26) 



(27) 
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Electric Wave Exx Propagation in Double Layered EM GL-Cloak 
In XY Plane Using 3D GL Modeling in GLGEO 
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E xx (r,r 8 ,t)\ r ^ Si ^0. (28) 

The electric intensity E xx (r, r S) t)\ reS ^ in (28) and 
E xx (r,r s ,t)\ reSi = 0. in (17) are an obvious contradic- 
tion. Therefore, we proved the Statement 2 that there 
exists no Maxwell EM wavefield can be excited by the 
nonzero local sources inside of the single layered cloaked 
concealment with normal materials. 



FIG. 1: (color online) At this moment of the time step 39dt, 
the most part of the front of the First electric wave, E xx ,i, 
propagates enter to the inner GL cloak layer, R\ < r < R2, a 
few part of the front of the E XXi i is located in right and top 
of the concealment; the front of Second electric wave, E X x,2, 
reaches the outer boundary r = R3 of the GL double layered 
cloak. 
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Electric Wave Exx Propagation in Double Layered EM GL-Cloak 




Horizontal X (Meter), Made In GLGEO 

FIG. 2: (color online) At time step 75dt, the 

Second electric wave, E XXt 2, is propagating inside of the outer 
layer cloak of the GL double layered cloak, R2 < r < R3, 
and around the sphere annular r — R2 and never pene- 
trate into inner domain, r < R2. It does disperse and split 
into the two phases around the sphere annular r = i?2- 
The First electric wave, E XXi i,is propagating inside of the 
inner sphere annular layer of the GL double layered cloak, 
Ri<r < R 2 . 



B. The EM Wave Excited By Point Source In The 
Concealment And Other Point Source In The Free 
Space Propagates Through The GL Double Layered 
Cloak 

The configuration of the GL double layered cloak mate- 
rial is described in the subsection A of this section. Two 
point sources are used to excite the EM wave propagation 
through the GL double layered cloak. The first point cur- 
rent source is located inside of the center sphere conceal- 
ment at (—0.12m, —0.12m, 0.0), by which the excited EM 
wave is named as First EM wave, its component E XXi i 
is named First electric wave, The second current point 
source is located in free space at (0.518m, 0.518m, 0.0) 
where is the right and top corner outside of the whole 
GL double layered cloak. The EM wave by the second 
source is named as Second EM wave. Its component 
E xx ,2, is named Second electric wave. The GL model- 
ing simulations of the EM wave excited by above two 
point sources propagation through the GL double lay- 
ered cloak are presented in the Figures 1-3. The two 
waves are propagating at time step 38<it that is shown 
in the Figure 1, at this moment, the most part of the 
front of the First electric wave, E xx> \, propagates en- 
ter to the inner GL cloak layer, i?i < r < R2, a few 
part of the front of the E XXt \ is located in right and top 
of the concealment. The front of Second electric wave, 
E xx .2, reaches the outer boundary r — R3 of the GL 
double layered cloak. In the Figure 2, at time step 75dt, 
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FIG. 3: (color online) At time step 98dt, one part of front the 
Second electric wave, E XX> 2, is propagating inside of the outer 
layer cloak , R2 < r < R3. It has around the sphere annular 
r = R2 and forward bending in the left of the sphere annular 
r = i?2- The First electric wave, E XX] i, is still propagating 
inside of the inner sphere annular layer of the GL double 
layered cloak, Ri < r < i?2- 



the Second electric wave, E xx ^, is propagating inside of 
the outer layer cloak of the GL double layered cloak, 
i?2 < r < i?3, and around the sphere annular r — R2 and 
never penetrate into the inner layer of GL cloak and the 
concealment , i.e. rJe.i?2- It does disperse and split into 
the two phases around the sphere annular r = R2 , the 
front phase speed exceeds the light speed; the back phase 
is slower than the light speed. The First electric wave, 
E XXt \, is propagating inside of the inner sphere annu- 
lar layer of the GL double layered cloak, R\ < r < R2. 
In the figure 3, at time step 98dt, one part of front the 
Second electric wave, E xx> 2, is propagating inside of the 
outer layer cloak , R2 < r < R3. It has around the 
sphere annular r — R% and forward bending in the left 
of the sphere annular r = R2 and never penetrate into 
the inner layer of GL cloak and the concealment, i.e. 
r < i?2 The First electric wave, E xxl , is still propagat- 
ing inside of the inner sphere annular layer of the GL 
double layered cloak, R\ < r < i?2- At time step I38dt, 
the Second electric wave, E xx $ has propagated outside 
of the whole GL double layered cloak, a small part of its 
wave front is located in the left and low corner of the 
plot frame which is shown in the Figure which is omit- 
ted, most part of front of the E xx ^ electric wave field has 
been out of the plot frame. The exterior EM wave out- 
side of the GL double layered cloak never been disturbed 
by the cloak and never penetrate enter the centre sphere 
concealment and inner layer of the GL cloak. At same 
time step, the First electric wave, E XXt i, is propagating 
inside of the the inner layer sphere annular of the GL 
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double cloak, R\ < r < R 2 . It can be very closed to 
the interface boundary r = R 2 , However, it can not be 
reached to the interface boundary r — R2 for any long 
time. 



Electric Wave Exx Propagation in Double Layered EM GL-Cloak 
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FIG. 4: (color online) At moment 21dt, one part of the front 
of the First electric wave, E xx ,i, propagates enter to the con- 
cealment; other part of the front is still propagating inside 
of inner GL cloak layer, Ri < r < i?2- One part of front 
of Second electric wave, E xx p, reaches the middle interface 
boundary r = R2 of the GL cloak; other part of the front has 
propagated in free space with disturbance. 



C. The EM Wave Excited By The Point Source In 

The Inner Layer Q,gli And Other Point Source In 
The Outer Layer Q.glo Propagates Through The GL 
Double Layered Cloak 

The 3D EM full wave excited by the point source 
in the inner layer and other point source in the outer 
layer of GL double Layered cloak are simulated by us- 
ing GL EM modeling. The simulations of the EM wave 
excited by the above sources through GL double lay- 
ered cloak are presented in the Figures 4. The config- 
uration of the GL double layered cloak material is de- 
scribed in the subsection A of this section. Two point 
sources are used to excite the EM wave propagation 
through the GL double layered cloak. The first point 
current source is located inside of the inner layer of the 
GL cloak at (—0.165m, —0.165m, 0.0), by which the ex- 
cited EM wave is named as First EM wave. Its com- 
ponent E XXj i is named First electric wave. The sec- 
ond current point source is located in outer layer GL 
cloak at (0.23m, 0.23m, 0.0). The EM wave by the sec- 
ond source is named as Second EM wave. Its compo- 
nent E XXj2 , is named Second electric wave. In the Fig- 
ure 4, at moment 21dt, one part of the front of the 



First electric wave, E xx ,\, propagates enter to the con- 
cealment; other part of the front is still propagating in- 
side of inner GL cloak layer, Ri < r < R 2 . One part of 
front of Second electric wave, E XXj2 , reaches the middle 
interface boundary r — R 2 of the GL cloak; other part of 
the front has propagated outside of the whole GL double 
layered cloak and in free space with disturbance. At the 
time step A8dt, the First electric wave, E xx \, has propa- 
gated through the concealment and whole front is inside 
of the inner layer of the GL cloak i?i < r < R 2 , its speed 
is smaller than the light speed. The Second electric wave 
never propagates into the inner layer of GL cloak, The 
part of the front of Second electric wave, E xx ^ 2 , is in- 
side of outer layer of the GL double layered cloak, R 2 < 
r < i?3, and being forward bending with speed lager 
than the light speed. Other part of the front has been 
propagating in free space with disturbance. The EM 
wave propagation image is presented in an omitted Fig- 
ure. In the other omitted Figure, at time step 68dt, the 
Second electric wave, E xx ^ 2 , is propagating in free space 
and outside of the whole GL double layered cloak, and 
never penetrate into inner domain, r < R 2 , i.e. never 
penetrate into the inner layer of GL cloak and the con- 
cealment. The First electric wave, E XXi ±, is still prop- 
agating inside of the inner sphere annular layer of the 
GL double layered cloak, R\ < r < R 2 . At time step 
98dt, the Second electric wave, E XXi2 , has been propa- 
gating out of whole GL cloak and out of the plot frame. 
Very small part of its front is located in left and low 
corner of the plot frame, the fihure is omitted. The 
First electric wave, E XXi \, is still propagating inside of 
the inner sphere annular layer of the GL double layered 
cloak, R\ < r < R 2 . It can be very closed to the interface 
boundary r = R 2 , However, it can not be reached to the 
interface boundary r = R 2 for any long time. 



VII. ADVANTAGES 

A. The EM GL Double Layered Cloak Is Robust 
For Invisibility 

The figure 4 clearly show that the wave front of the Sec- 
ond electric wave has propagated outside of the cloak and 
go to free space with disturbance. The results reminders 
us that if only single outer layer cloak floLO is adopted, 
and there is a little crack loss on the inner boundary sur- 
face dflcLO-: some EM or current source inside of the 
^glo w iU excite the EM wave propagation go out to free 
space and expose the cloak immediately. The GL double 
layered cloak overcomes the weakness that is also shown 
in the figure 4. The wave front of the First electric wave, 
which is excited by a point source inside of the inner layer 
cloak Qqli , is always propagating inside of inner layer 
cloak Qgli or concealment fl C oni and never propagate 
outside of the interface annular Si. Therefore, the EM 
GL double layered cloak is robust for invisibility. 
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The EM GL Double Layered Cloak Is Complete 
Invisible 



The figure 1-3 clearly show that the outer layer cloak of 
the GL double layered cloak has the invisibility function, 
while the inner layer cloak has fully absorption function 
to absorb the EM wave excited from local sources inside 
of the concealment. When the source is located outside 
of the GL double layered cloak, the excited EM wave field 
propagation outside of the double layered cloak is as same 
as in free space and never be disturbed by the cloak; also, 
the exterior EM wave can not penetrate into the inner 
layer and concealment of the GL double layered cloak. 
When local sources are located inside of the GL double 
cloaked concealment with the normal EM materials, the 
excited EM wave is propagating under Maxwell equation 
governing, it is complete and successively absorbed by the 
inner layer cloak of GL double cloak and never propagate 
to outside of the inner layer of the GL cloak, moreover, 
the EM wavefield in concealment never be disturbed by 
the cloak. 

Using the GL method theoretical analysis, the state- 
ment 2 in section 4 is rigorous proved. It states that 
"there exists no Maxwell electromagnetic wavefield can 
be excited by nonzero local sources inside of the single 
layer cloaked concealment with the normal EM mate- 
rials". The invisibility of the single layered cloak and 
existence of Maxwell EM wave field excited by the local 
sources inside its concealment is inconsistent. Provide 
only single outer layered cloak is adopted. The EM field 
excited by local sources inside of its concealment with 
normal materials does not satisfy the Maxwell equation. 
The EM chaos phenomena, which is divorced from the 
Maxwell equation governing, may damage devices and 
human inside the concealment, or may degrade the in- 
visibility of the cloak. Invisibility function of the single 
layered cloak is not complete. The EM GL double lay- 
ered cloak overcomes the drawback of the single layered 
cloak. The GL double layered cloak has the complete 
sufficient invisibility function. 



C. The EM GL Double Layered Cloak Can Be 
Double Ellipsoid Annular 

The EM GL double layered cloak can be extended to 
have double ellipsoid annular and other geometrical dou- 
ble layered closed strips. 



D. Frequency Band 

Many simulations and theoretical analysis by the GL 
method show that the idea EM GL double layered cloak 
is of the invisibility function for all frequencies. However, 
the practical material has some loss. The frequency band 
will be depended on the rate of the material loss. Because 
the EM GL double layered cloak is robust and complete 



cloak, it has three radius R2, and R3 can be chosen, 
and GL method reduced numerical frequency limitation 
in FEM and FD scheme, a reasonable wide frequency 
band of the GL double layered cloak for low loss rate will 
be optimization and reported in next paper. 



E. Advantages Of The GL Method 

The GL EM modeling is fully different from FEM and 
FD and Born approximation methods and overcome their 
difficulties. There is no big matrix equation to solve in 
GL method. Moreover, it does not need artificial bound- 
ary and absorption condition to truncate the infinite do- 
main. Born Approximation is a conventional method in 
the quantum mechanics and solid physics, however, it is 
one iteration only in whole domain which is not accu- 
rate for high frequency and for high contrast materials. 
The GL method divides the domain as a set of small 
sub domains or sub lattices. The Global field is updated 
by the local field from the interaction between the global 
field and local subdomain materials successively. Once all 
subdomain materials are scattered, the GL field solution 
is obtained which is much more accurate than the Born 
approximation. GL method is suitable for all frequency 
and high contrast materials. 

Moreover, the GL method can be meshless, including 
arbitrary geometry subdomains, such as rectangle, cylin- 
drical and spherical coordinate mixed coupled together. 
It is full parallel algorithm. These advantages of the GL 
method to overcome historical difficulties have been de- 
tailed described in the paper [1]. The theoretical founda- 
tion of the GL method is described in the paper [2] The 
GL EM method consistent combines the analytical and 
numerical approaches together and reduced the numeri- 
cal dispersion and numerical frequency limitation. The 
GL method has double abilities of the theoretical anal- 
ysis and numerical simulations that has been shown in 
this paper. 

The 3D GL simulations of the EM wave field through 
the single and multiple sphere, cylinder, ellipsoid, and 
arbitrary geometry cloaks in single layered and double 
layered profile show that the GLT and GLF EM modeling 
are accurate, stable and fast. It saves much more storages 
than the conventional methods. In general, only 10 to 50 
minute are needed to run the 3D EM wave field through 
the cloaks with 64 to 128 frequencies in the PC. The 
high performance GL parallel algorithm in PC cluster 
and super parallel computer is very fast and powerful to 
simulate complex and large scale physical and chemical 
process. 

A double layer cloth phenomenon to prevent the 
GILD inversion [5-7] detection has been observed 
in paper [8] in 2001 which is published in SEG online 



http://www.segdl.org/journals/doc/SEGLIB-home/dci/searchDCL. 
After the event, we effort improve GILD [5-7J and devel- 
oped a novel and effective Global and Local field (GL) 
modeling and inversion to study the meta materials, 
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periodic photonic crystals and condense physics etc. 
wide physical sciences. 3D GL EM modeling and 
inversion [9] and computational mirage [10] have been 
presented in PIERS 2005 and published in proceeding of 
PIERS 2005 in Hangzhou, which can be downloaded from 
http://piers.mit.edu/piersproceedings/piers2k5Proc.php. 
We developed 3D FEM for the elastic mechanics first 
in China in 1972 which has been cited and recorded in 
[17]. Our 3D FEM paper has been published in [18] in 
Chinese. We are deeply and clearly to know the merits 
and drawbacks of FEM and its serious limitation and 
difficulties to simulate high frequency wave propagation 
through dispersive materials. The GL method overcome 
the drawbacks of FEM and FD methods. The history 
of development of our 3D FEM [18], novel inversion [7], 
GILD [5] and GL method [1-3] has been described in 
[2]. The 3D and 2D GL parallel software is made and 
patented by GLGEO. The GL modeling and its inversion 
[1-3] and GL EM quantum field modeling are suitable 
to solve quantization scattering problem of the electro- 
magnetic field in the dispersive and loss metamaterials, 
cloaks and more wide anisotropic condense materials. 



VIII. CONCLUSIONS 

Many simulations of the EM wave propagation through 
the GL double layered cloak by the GL modeling and the- 
oretical analysis verify that the EM GL doubled cloak is 
robust cloak and has complete and sufficient invisibility 



functions. Its concealment is the normal electromagnetic 
environment. The outer layer cloak of the GL double 
layered cloak has the invisible function, the inner layer 
cloak has fully absorption function. When the source is 
located outside of the GL double layered cloak, the ex- 
cited EM wave field propagation outside of the double 
layered cloak is as same as in free space and never be dis- 
turbed by the cloak; also, the exterior EM wave can not 
penetrate into the inner layer and concealment. When 
sources are located inside of the GL double cloaked con- 
cealment with the normal EM materials, the excited EM 
wave is propagating under Maxwell equation governing, 
it is complete absorbed by the inner layer cloak of GL 
double cloak and never propagate to outside of the in- 
ner layer of the GL cloak, moreover, the EM wavefield in 
concealment never be disturbed by the cloak. The EM 
GL double layered cloak has advantages to overcome the 
drawback and difficulty of the single layered cloak. 

The GL method is an effective physical simulation 
method. It has double abilities of the theoretical analysis 
and numerical simulations to study the cloak metamate- 
rials and wide material and Field scattering in physical 
sciences. 
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